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Abstract 

We write down an action for a charged, massive spin two field in a fixed Einstein back- 
ground. Despite some technical problems, we argue that in an effective field theory framework 
and in the context of the AdS/CFT correspondence, this action can be used to study the 
properties of a superfluid phase transition with a ci-wave order parameter in a dual strongly 
interacting field theory. We investigate the phase diagram and the charge conductivity of 
the superfluid phase. We also explain how possible couplings between the spin two field and 
bulk fermions affect the fermion spectral function. 
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1 Introduction and summary 

In the context of gauge gravity duality [IJ 12} E], a holographic superconductor has come 
to mean a strongly interacting field theory that undergoes a superconducting or superfluid 
phase transition and that has in addition a dual gravitational description. See jH [6] 
for reviews. Holographic superconductors are interesting because they describe condensed 
phases of strongly coupled, planar, gauge theories. Also, they may give some insight into 
real world strongly interacting superconductors and superfluids such as the cuprates and 



helium-3. To date, most of the holographic examples have had either s-wave or p-wave 
order parameters (three exceptions are [7J El E]). Because of a possible relationship with the 
cuprates, in this paper we explore a holographic toy-model for the condensation of a <i-wave 
order parameter. 

A toy-model for an s-wave holographic superconductor was first constructed in [TO] fol- 
lowing the observation in [TT] that in AdS space black holes can have scalar hair. A more 
realistic construction of an s-wave holographic superconductor dual to an orbifold of M = 4 
SYM was discussed in |12j . See [13] for a similar construction using an M-theory truncation 
and [H] for a p-wave construction using probe D-branes. 

Comparing a holographic superconductor to a high T c superconductor or other materials 
is tempting but difficult: For a gauge theory such as M = 4 SU(N) Super Yang-Mills, the 
dual gravitational description is most useful in its strongly coupled and large N limit. Moving 
away from this limit is challenging. The role of supersymmetry and the absence of a lattice 
are two other important issues that must be confronted before any detailed comparison can 
be made. Nevertheless, the need for non-perturbative insight in studying high T c materials 
[15] suggests an exploration of holographic methods. 

In the holographic setup, a G?-wave order parameter is dual to a charged massive spin two 
field propagating in an asymptotically AdS geometry. While the action for a neutral, massive 
spin two field in a flat background has been known for over seventy years [16J, finding an 
action for a charged spin-two field in a curved spacetime is an unsolved problem. In section 
[2] we discuss the perils and pitfalls in constructing such an action. We settle for a relatively 
simple action with some known problems that we believe can be fixed through the addition 
of higher dimensional operators. Thus, we work in an effective field theory framework. 

With the action for a charged massive spin two field in AdS space in hand, we can solve 
the bulk equations of motion and study the resulting condensed phase in the boundary 
theory. This is carried out in section [3j For a particular ansatz in which only a single, real 
component of the spin two field is turned on, and after an appropriate field redefinition, 
we observe that the equations of motion for the spin two field are identical to those of the 
charged massive scalar field studied in [TO]. Thus, we are guaranteed that when turning on 
a chemical potential there will exist a critical temperature below which the spin two field 
condenses. 

We proceed to compute the conductivity of the d-w&ve superconductor in section [4] and 
to study fermion correlators in this background in section [5} As was the case for the massive 
charged scalar, we find that the conductivity associated with the superconductor develops a 



gap-like feature (depicted in figure [2j) but it also has a sharp (^-function spike indicative of 
a bound state. In studying the spectral function for fermions, we find that for a particular 
coupling between the fermions and the spin two field, the spectral function has a band-like 
structure in addition to the d-wave gap and Fermi-arcs observed in [H]|j 

2 Spin two fields 

One way of generating a <i-wave condensate in the boundary theory is to have a massive, 
charged spin two field condense in an asymptotically AdS^+i geometry. By a massive spin 
two field in d + 1 spacetime dimensions, we mean a field that transforms locally in the 
(d + 2)(d — l)/2 dimensional irreducible representation of the little group SO(d) of the 
Lorentz group SO(l, d). If such a field is represented by a symmetric tensor (p^ then there 
should be enough constraint equations to eliminate d + 2 of its (d +l)(d + 2)/2 components. 
In flat space, the Fierz-Pauli Lagrangian [16] does this job. It is given by 



H 



- d^^ifT + 2<p^ - 2^d^ + d^ip - m 2 {y^ v - ^ 



where we have introduced the notation <p p = d^ip^p and if = tp't. The equations of motion 
(EOM's) and constraints following from (flj are 

= (□ -m 2 ) w , (2a) 

= daifT , 

(2b) 
which indeed give the correct number of degrees of freedom. We reproduce pi) in appendix 

EH 

Naively, one might think that covariantizing (fTl) and introducing a minimal coupling to 
a U(l) field will yield a consistent action for a massive charged spin two field in a curved 
background. Unfortunately such a Lagrangian does not produce appropriate constraint 



equations as in (2b) and the spurious propagating degrees of freedom give rise to pathologies 
such as ghosts, loss of hyperbolicity and faster than light travel. There have been several 
attempts to deal with these issues in the literature. We mention two. The Buchbinder- 
Gitman-Pershin (BGP) Lagrangian [19] describes a neutral massive spin two field in an 
Einstein manifold, and has the appropriate number of degrees of freedom that propagate 



L See [T5] for a related discussion of Fermi arcs in p-wave superconductors. 



causally. The Federbush Lagrangian constructed in [20J describes a charged massive spin 
two particle propagating in flat space. It has the correct number of propagating degrees of 
freedom, but also has the unwanted feature of generating superluminal modes [21] . 
The Lagrangian density we will use has the following form: 



C = -|ZV/Vl 2 + 2|IWT + \ D M 2 ~ [D^ u D uV + c.c] - m 2 (|¥v| 2 - \y\ 2 ) 



1 „, ,o . „ .... ,. 1 



+ 2R iiVpX ^ H T x - R^Vx - J^ R \V\ 2 - iqF^f^Vx ~ j V" 



(3) 



where we have introduced D p = V M — iqA p when acting on (p pu and now ip p = D^ip^. The 
EOM's which follow from (J3J) are 

= (□ -m 2 )ip pv - 2D (p Lp u) + D {il D v) (p - g pu [(Q -m 2 )(p - D p ip p ] 

+ 2R ppuX ^ x - 9^-^P - i\ {F w <f>$ + F vp <pfi (4) 



d p f^ = r 



where 



r = vplpiiy^ - D"<ft*>) + i(& - D a( p*)(^ - gT<p) + h.c. . (5) 

The Lagrangian (|3J) has the following desirable properties: 



It has familiar limits. For a neutral, non-interacting spin two field in flat spacetime, 
it reduces to the Fierz-Pauli Lagrangian [IE] . For a neutral spin two field in a fixed 
Einstein background, C reduces to that of BGP [19]. For a charged spin two field in 
flat spacetime, C reduces to the Lagrangian of Federbush 



• It is ghost-free, and for generic values of F pu , it describes (d + 2)(d— l)/2 propagating 
degrees of freedom. We demonstrate this fact in Appendix [A] 

• The Lagrangian is unique in an appropriate sense: Consider a spin two Lagrangian with 
all possible operators up to dimension d + 1 and quadratic in tp pu . If we choose the 
couplings of this Lagrangian to be different from those in (pi), then some of the would-be 
constraint equations become dynamical, introducing ghosts. We provide more details 
about the uniqueness properties of ^ in appendix |AJ 

At the same time, our Lagrangian has a number of limitations: 

• As was the case for the BGP Lagrangian, we are restricted to work in a fixed back- 



ground spacetime that satisfies the Einstein condition 

R-nv = j—^ g^u • (6) 

Otherwise, the constraint equations will not be satisfied and we will have too many 
propagating degrees of freedom. In the context of holographic superconductors, this 
restriction forces us to work in the probe limit [10J where the spin two field and gauge 
field do not backreact on the metric; since the metric will not be perturbed by the 
matter content of the theory, it will automatically satisfy (|6]) . 



Similar to the EOM's following from the Federbush Lagrangian [20], the EOM's follow- 
ing from (J3J), for generic values of F^ u , are either non-hyperbolic or lead to non-causal 
propagation [21J. Fortunately, these effects are small in the sense that they disappear 
in the limit where the quantities 

9 ' 5 \' ) 

vanish. (Underlined indices are in a local vielbein basis.) As we summarize in section 



3.3, there are a number of reasons to believe that one can correct these problems by 



adding to the Lagrangian terms that are higher order in -Fi^rl Such a series expansion 



is sensible only in a frame where these violations of causality are small. In section 3.3 



we will discuss in what regime of parameter- space these corrections can be made small. 

3 The bulk dual of a d-wave condensate 

In this section we adapt the standard AdS/CFT map to a massive spin two field, and study 
the gravitational solution corresponding to a holographic (i-wave superconductor. 

3.1 The conformal dimension of the boundary operator. 

Since AdS space is an Einstein manifold, we can use the Lagrangian C defined in (J3J) to 
describe a charged spin two field propagating in AdS space. The gauge gravity duality 



2 One of the reasons which leads us to believe that these corrections exist is the existence of the Argyres- 
Nappi Lagrangian 22J. The Argyres-Nappi Lagrangian is a consistent causal Lagrangian for a massive 
charged spin two particle in a constant electromagnetic background and 26 flat spacetime dimensions. Further 
indication that higher order terms containing F^ are needed comes from the fact that the Lagrangian pi, 
when rewritten in terms of Stiickelberg fields and in a gauge that makes the kinetic terms canonical, contains 



such terms [23j|24]. We will come back to this point in section 3.3 



implies that this Lagrangian also describes a spin two operator O mn on the boundary theory. 
Since £ was not derived from a higher dimensional string action, it is not clear what the 
boundary theory dual to (pi) is, or if it exists as a proper field theory. In this sense the action 
(|3| should be thought of as a toy-model. Assuming the boundary theory exists, we can use 
the AdS/CFT dictionary described in [21 E] to make some general observations regarding the 
dual boundary theory spin two operator. The following analysis of the correspondence for 
a spin two field generalizes the well known duality between the energy momentum tensor of 
the boundary theory and the graviton — a neutral, massless spin two field in the bulk. 

Consider the spin two Lagrangian ^ for a neutral spin two field propagating in an AdS 
background with line element 

L 2 

ds 2 = — ( - dt 2 + dx 2 + dz 2 ) . (8) 

As written, it is not clear from (|3| how to define the mass of the spin two field, in the 
sense that one can absorb couplings between the spin two field and the curvature into a 
redefinition of m 2 [25] . We resolve this ambiguity by defining the mass squared of the spin 
two field, m 2 , such that when m? = the action has an enhanced gauge symmetry of the 
form ip^ v — > ip^ v + 5(p^ u with 5ip^ y = D^ u + D u ^, where ^ is an infinitesimal parameter. 
This gauge invariance is unique to the graviton, and is associated with diffeomorphisms. 

In what follows we will be setting m 2 =£ 0. When m 2 ^ and m 2 ^ dtd+D ^i ^ ne equations 
of motion and constraints derived from C are similar to the flat spacetime case pi): 

= ( □ -m 2 ) w + 2R lxXvp <p x e (9a) 

(9b) 



o- '«" 



Curiously, when m 2 = ^Ay-R the spin two action has an enhanced gauge symmetry of the 
form 5(ffj, u = D^D U £ + g^ d i^i x \ £■ This is related to the existence of a partially massless 
spin-2 field |2B]. 

For operators with spin s > 0, bulk fields and boundary operators have a different number 
of components. In what follows, we will argue that the equations of motion for a bulk spin 
two field (f^u have two linearly independent solutions, and that the number of integration 
constants is twice the number of components of a spin two operator O mn in the boundary 



theory. In 3.2 we will identify these solutions with the source and expectation value of the 



spin two operator. 



The explicit form of the constraint equations (9b) is 



[(d - 1) - zd z ]ip mz = zr] kn d k Lf nm 
[(d- 1) - zd z ](p zz = zr] kn d k (p nz (10) 

where roman indices run over the boundary directions 0, . . . , d — 1 and rj kn is the Minkowski 
metric. We observe that the components of ip ul/ with one or more legs along the radial 
direction are completely determined by the components of ip^ with both legs along the 
boundary directions. Furthermore, the near boundary series expansion of (p^ z starts off with 
a higher power of z than the near boundary series expansion of ip m n- Therefore, the last 



equation in (10) implies that not only is ip UL/ traceless, but that the leading term in a near 
boundary series expansion of p> mn is also traceless. Given the relation between (O mn ) and 
tp mn which we make precise below, tracelessness of the leading term of a near boundary series 
expansion of (f mn ensures that the boundary operator O mn is also tracelessjj 

The equations of motion for the remaining components of the spin two field are 

= -[2(d-2)+L 2 m 2 ]ip mn -2r] mn T] kp ip kp + z d - 3 d z (z 5 - d d z ip mn )+z 2 d 2 ip mn -Azd {m if n)z . (11) 



These equations have two solutions whose near boundary behavior is given by (p mn ~ z A 2 
and (p mn ~ z d ~ A ~ 2 , where A is defined through 

m 2 L 2 = A(A - d) . (12) 

The parameter A should be identified with the dimension of the operator O mn as we now 
explain. Using the AdS/CFT dictionary, the coefficient of the z A ~ 2 term can be associated 
with an expectation value, (O mn ), while the coefficient of the z d ~ A ~ 2 term serves as a source, 
Omn- Since (p nu has definite scaling behavior under z — ¥ z/X, the ratio (O mn )/Om n must 
have scaling dimension 2A — d. This ratio, through the theory of linear response, is the 
Fourier transform of a retarded Green's function which, for an operator of dimension A, will 
also have scaling dimension 2A — dr\ As expected, when A = d (which is the protected 



3 The equation of motion for the spin two field is a second order differential equation. Tracelessness of the 
leading term in a series expansion of the solution applies to both of the two linearly independent solutions. 
Thus, both the source term and the expectation value of the dual operator will be traceless. 

4 As a rule of thumb, the two near-boundary behaviors of the boundary components $ mi ... m of a generic 
tensor field V^...^ with covariant action and in the standard coordinates (|8J), are z A ~ p and z d ~ A ~ p , where 
A is the dimension of the boundary operator C mi ... m . This scaling can be motivated by noticing that the 



dimension of the boundary theory stress tensor) then m 2 = 0. 

The dimension of a spin two primary operator in d spacetime dimensions is constrained by 
unitarity and conformal symmetry. An analysis of the positive-energy unitary representations 
of the three dimensional conformal group 5*0(3, 2) was first carried out in [27]. The analysis 
has been generalized to the four dimensional conformal group 5*0(4,2) in [2E] and to d 
dimensions in [29]. In short, the analysis in [271 EE1 EH] states that the positive-energy 
representations of SO(d,2) are classified by their quantum numbers in 50(2) x SO(d), 
which are the dimension A of the primary operator in the module and its spin. Unitarity 
sets a spin-dependent lower bound on A. In particular, a totally symmetric and traceless 
tensor operator O mi ,,, ms of spin s > 1 must satisfy A > d — 2 + s. When this bound is 
saturated, then O mi ,„ mg is a conserved current in the sense that d mi O mi „, mg = 0. For spin 
two operators the unitarity bound implies that A > d. Thus, our bulk spin two fields should 
have a positive mass squared, m 2 > 0. In (2H1 EHl EI] (see also [32]) it was shown that also 
from the bulk point of view, stability implies that m 2 > 0. In appendix [B] we discuss how to 
obtain such a bound on the mass by generalizing the analysis of Breitenlohner and Freedman 
[33]. 

3.2 Generating a d-wave condensate 

The spin two action we consider describes the correct number of degrees of freedom for 
a massive spin two field only on an Einstein background. Therefore, we need a control 
parameter such that the bulk stress tensor can be made negligible when compared to the 
cosmological constant in the Einstein equations. This regime of the theory, called the probe 
limit in [10J, is achieved by taking the charge of the spin two field q to be large. Indeed, 
rescaling the spin two field and the gauge field by a factor of q, (p^ = (p^vjq and A^ = A^/q, 
the Lagrangian (J3J) gets rescaled by a factor of q 2 , £ = C/q 2 (where C is independent of q). 
In the limit in which q is infinite but A^ and tp^ are kept fixed, the bulk stress tensor goes 
to zero. We will always work in this limit, expressing boundary quantities in units of qfi, 
where /x is the chemical potential in the boundary theory. 

To look for a thermal state on the boundary theory in which the spin two field presumably 
condenses, we need to solve the equations of motion Q in a black hole background which 



boundary and bulk metrics are related by the rescaling g hnd = z 2 g hulk \, ,. Thus, the identification of tensor 
quantities expressed in a boundary or bulk frame requires them to be rescaled accordingly. 



satisfies the Einstein equation (J6J). One such geometry is given by the black brane solution 

ds 2 = ^ ( - f(z) dt 2 + dx 2 d _ x + ^-A (13) 

where 

f(z) = l-(-Y. (14) 

\z h / 

The black hole horizon is located at z = Zh, while the conformal boundary of the spacetime 

is located at z — 0. The Schwarzschild temperature of this black hole is 

T =T~- ( 15 ) 

We consider an ansatz where ip^ v and A^ depend only on the radial coordinate z and 
where only the space components of ip^ are turned on. For d = 3, since (p^ u is a complex 
field, we can use a spatial (xy) rotation combined with a U(l) gauge transformation to 
choose one of the two components of ip, ip xx = —<p yy or ip xy , to be real and the other one to 
be imaginary. However, similar to what was found in [10], it is consistent to turn on a single 
real component of ip and also to set all the components of the gauge field except for A t to 
be zero. Our ansatz is then 



L 



2 



A = A^dx^ = <j)(z)dt , (p xy (z) = —^ip(z) , (16) 

with all other components of (p^ u set to zero, and and ip real. The more general ansatz in 
d = 3 with one mode real and the other one imaginary could lead to interesting physics and 



we leave it for future work. The ansatz (16) satisfies <p = ip^ = F^pip? = 0. Instead of turning 
on <p xy in (16) we could have considered a non- vanishing value for <p xx - yy = <p xx = —<p yy - 
These two ansatze are equivalent under a 7r/4 rotation. Indeed, under a counterclockwise 
rotation in the xy-plane by an angle 6, the vector ( Vx ^~J y ) transforms under the action of 

the matrix 

f cos 29 -sin2fl\ , x 

Re=[ . 17 

\sin2# cos 26 J 

Note that our ansatz continues to satisfy the equations of motion when d > 3, but the most 
general ansatz will require turning on more components of ip^. 



Using the ansatz (16), the EOM's for (ft and ip are 







Ul I ^ _ ^ xl 1-^2 



+ 



z 2 f 



2x1 



"^l.f-^W? 



2 r2 



m L 



1ft 



(18) 
(19) 



where a prime denotes a derivative with respect to the radial coordinate z. Near the asymp- 
totically AdS boundary the series expansion for ip takes the formjj 



ij,( z ) = z d ~ A [ij {s) + 0(z)] +z A 



l2A-d 



+ 0{z) 



(20) 



where A is the conformal dimension of the boundary spin two operator defined in (12). 



The O(z) terms represent sub leading corrections in powers of z which can be computed 
perturbatively from the EOM's. The values ip^ and (O) are integration constants which 
need to be determined. 

The boundary conditions we impose on the spin two field are that the modulus (p f _ iu tp fiu is 
finite at the black hole horizon, which implies that ip is finite. We also require that ip^ = 
so that the spin two field is not sourced. The gauge field A t is required to vanish at the 
horizon. We interpret the value of A t at the asymptotically AdS boundary as the chemical 



potential /x. The equations (18) and (19) are identical to the EOM's for the Abelian Higgs 



model in AdS^+i [ID]. Thus, a condensed solution exists below a certain critical temperature 
T c : besides the "normal phase" solution 



A 



(z)dt 



fJ 



Zh 



d-2- 



(It 



^ = 



(21) 



which exists for all values of the chemical potential, an additional solution with ip ^ 
exists whenever T/q/i is small enough. If we consider the chemical potential as fixed, then 
the critical value of T/qfj, defines a critical temperature T c below which the spin two field 
condenses. 



To determine T c , we look for the temperature at which the normal phase solution (21) 



becomes unstable to perturbations in the sense that ip develops a normal mode. There exist 
several strategies for determining T c . For large m 2 the problem becomes classical, and the 
critical temperature can be approximated by determining the condition for a global minimum 
away from the horizon in the potential of a charged, massive particle in the normal phase. 



3 Up to possible logarithms for special values of A. 
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We examine this method in Appendix [C] To determine T c for other values of the mass, and 
to find an explicit solution for T < T c , we must resort to numerical techniques. 



With a numerical solution to (18) and (19) in hand, one may compute the expectation 



value of the boundary operator O mn dual to the spin two field (p mn via the standard AdS/CFT 
rules with flat Cartesian coordinates on the boundary [21 E] : 

(O mn ) = lim— -r-, (22) 

where (fm n = lim 2 _>o z A ~ d+2 (f mn (z) is the source term, and S is the on-shell action plus pos- 
sible boundary terms£j Instead of working out the boundary counter terms for the spin two 
field and computing (O mn ) directly, we will use the similarity between the spin two system 
and the Abelian Higgs model in AdSd + i which was mentioned earlier. It is straightforward to 
show that the on-shell Abelian Higgs action discussed in [10] is identical to the one following 
from PJ and the ansatz (16). The only counter term which is required to render the on-shell 



Abelian Higgs action finite is a mass counter term. One can show that such a counter term 
is equivalent to adding a boundary counter term proportional to |</? M *,| 2 to ^. A typical 
dependence of (O xy ) on the temperature, obtained numerically, is shown in figure [I] 

We conclude by stressing that in the probe limit it is not possible to reach zero tempera- 
ture when m 2 > 0. In the probe limit, the background generated by the spin two field is the 
same as that of the Abelian Higgs model. The zero temperature limit of the Abelian Higgs 
model with backreaction was analyzed in (36J [37] where it was shown that when m 2 > the 
backreacted geometry has Lifshitz symmetry. In the probe limit, in contrast, the geometry 
has conformal symmetry when the temperature is zero. Thus, at zero temperature the back- 
reaction of the fields produces a qualitative difference on the geometry. Still, for any positive 
temperature T, it is always possible to take q large enough such that the stress tensor is 
negligible]^] 

3.3 Regime of validity of the effective Lagrangian 

In equation ([7]) we claimed that the field strength and its derivatives must be small when 
measured in units of the mass of the spin two field in order for our spin two action pi) to be 



6 In [34] a different normalization for ifj was suggested. This will result in a rescaling of (O xy ) by an 
overall A-dependent multiplicative factor. 

7 An additional problem with the T — > limit at non-zero density is that, although the action near T c 
has some claim to universality, at lower temperatures irrelevant operators become important and then the 
simple action we have written down needs to be corrected. 
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Figure 1: (Color online) The order parameter (O xy ) in d — 3 spacetime dimensions as a 
function of the temperature for various values of A which decrease from top to bottom. The 
chemical potential \x and the charge density p were obtained numerically using standard 
methods. The units we use to measure (O xy ) are similar to those used in [35] . 



valid as an effective theory. More precisely, we require all the components of F\ v in a local 
vielbein basis (which we denote F\^} to be smalljj In the AdS-Schwarzschild background 
where A t is given by (12X1) the only non-vanishing component of the field strength is 



qFu (d - 2)qpz 



d-i 



m z 



(23) 



A(A - d)z{ 

At any non-zero temperature the coordinate z can be no larger than Zh, yielding the condition 

(d - 2)q/iz h 



A(A-d) 



<1. 



(24) 



For large m 2 and d = 3 appendix [C] provides an estimate of the critical temperature below 
which the spin two field condenses. In order for the phase transition to occur in a regime 



where we can still trust the spin two action we find that (24) and (87) imply 



A>1 



(25) 



8 This statement is frame-dependent. The existence of a frame where all field strength components are 
small assures that the difference between our effective action and the full causal action (that we assume to 
exist) organizes in a series expansion. 
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By the same analysis, terms involving the gradient of F^ u in (l7| will be suppressed by an 
additional power of A. 

In principle effects of non-hyperbolicity and non-causality can be made small enough by 
working with a spin two field of very large conformal dimension. As we now explain, in 
practice it is difficult to access the solutions numerically when A is large. The requirement 
that the boundary spin two operator is not sourced implies that the leading z d ~ A ~ 2 term in 
the series expansion of the bulk spin two field must vanish. The expectation value of the 
condensate is proportional to the z A ~ 2 term in the series expansion of the bulk spin two field. 
Thus, when shooting from the horizon to the boundary, one needs to set the source term to 
zero to very high precision in order for the value of the expectation value to be accurate. 
In our numerical work, we chose a moderate value for A where we feel confident that the 
numerical error is small. We hope that the results will not qualitatively change with a larger 
A or with corrections to the action. 

Unfortunately, our expansion in inverse powers of A implies that dimension six operators 
of the schematic form 

— 2 Dip* Dip (26) 

will not be suppressed relative to the dimension four operator Fip*ip that we did consider 
in constructing pi). One way to motivate the addition of such dimension six operators is 
that they appear in the Argyres-Nappi action [22]. Using bosonic string theory, Argyres 
and Nappi found a consistent action for a spin two field in flat spacetime and constant 
electromagnetic background. If one expands the Argyres-Nappi action to linear order in 
F^ u , one finds terms precisely of the form p6J)r] 

The terms in (26) have a benefit. Using the effective Lagrangian pi), the gyromagnetic 



coupling g of the spin two field takes a surprising value of 1/2. Physical arguments suggest 
that g = 2 [3B1 EH]. The reason for this mismatch is that the gyromagnetic coupling is 
affected by irrelevant corrections and cannot be read off from our leading order Lagrangian; 



corrections of the form (26) alter g and the value g = 2 can be easily obtained. 

So why don't we include these dimension six operators? The honest answer is that we 
currently do not have a prescription for fixing their coefficients uniquely and that they make 
the calculations much more complicated. However, we can show that, for the most part, 



our results will be independent of these terms. With the ansatz (16), we have the relations 



9 The action is only consistent in d = 26 dimensions. For d ^ 26, a constraint equation becomes dynamical 
through the appearance of terms quadratic in F^. Thus, the action has the correct number of degrees of 
freedom and is causal only up to 0(F 2 ). It is not presently known whether a modification of the original 
action can work in arbitrary dimension. 
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F^v^p — and ip^ = 0. These equalities imply that the problematic dimension six operators 



in (26) won't affect the condensed solution described in section 3.2 Similarly, these operators 
won't affect computations of fermion correlators in the condensed phase performed in [UJ 
and discussed in more detail in section [5j However, for the conductivity calculation below, 
these operators will probably play a role. Given the added complexity of the equations one 
must solve, we content ourselves by studying the simpler action pi) and leave a detailed 



investigation of the operators (26) for the future. 



4 Conductivity 

The optical conductivity cr(uj) associated with the spin two condensate is a measure of the 
response of the (weakly gauged) U(l) current to a time varying external source, 

J = a-E, (27) 

where E = E e~ tu)t with E a constant electric field and a a tensor. To compute the 
conductivity in a holographic framework we turn on a source for the current J dual to the 
gauge field in the bulk. Using the holographic dictionary [2j |3], when the metric on the 
boundary theory is Minkowski, this current is given by 

(J<) = ^4 = (d - 2) At 2) , (28) 

where S is the on-shell action including possible counter terms and the A\ are given by a 
near boundary series expansion of the gauge field A: A4 = ^™ A- 2". The indices i, j, k 
run from 1 to d — 1. The leading coefficient A\ plays the role of the source term. 

The bulk equations of motion couple linear fluctuations of the complexified gauge field 
Ai to other complexified vector fields. In the probe limit these are tp ti and y? ti . This coupling 
between the vector modes implies that turning on an electric field will source not only the 
current J but also the fluxes K t = (Ou) and K* = (0*,). More explicitly, 

K = r-E (29) 

with t a rank two tensor, and there is a similar equation relating K* to E. Moreover, the 
source terms for (Ou) and (O u ) will source not only the fluxes K and K* but also the current 
J . Thus, in addition to the conductivity and the two aforementioned {d— 1) X {d — 1) tensors, 
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we have six additional (d — 1) X (d — 1) tensors describing the response of K, K* and J to 
the source terms for (On) and (O u ). 

In what follows we will focus on d = 3. To compute the frequency dependent conductivity 
in the x-direction, we consider linear, time dependent fluctuations of the complexified fields 
A x , (f ty , (p* y , ip zy and (p* zy . At linear order, the coupled algebro-differential equations for the 
e -«wt com p 0ri ent of these fluctuations are 

= K + J K + J 2 A* + 0- 2 [(<" - W)<Pt V ~ (« + mVtu] 

- '-f « - <flfj + j- f W ~ *f) {< ~ ^v) > ( 30a ) 

// 2 . 2f + m 2 L 2 2 qu + 2q 2 (p % r -, 

= Ptj, + - Vtw ^27 ^ + L — 7^7 — ^ + 2 L ^ + WWtv + g< ^ ^*w ' 



(30b) 



= [{u + g0) V - m 2 L 2 /] ^ + j tfg.ftK, + ^ WA 



- i(u + q4>)z 2 ip' ty - - [4(w + g0)^ + g^ 2 ] ^ . (30c) 

In all but the first equation we expressed /' and /" in terms of /. Notice that since tpty, 
(pl y , (p zy , (p* zy have to be treated as independent modes, the equations for the starred modes 
are not the complex conjugate of the unstarred ones — rather they are obtained by complex 
conjugation and the additional transformation u> — > —u. The functions (p zy and (p* can 



be eliminated from the first two equations using (30c), leaving three coupled differential 
equations for A x , ip ty and ip* ty . Since A y does not couple to this set of fluctuations, we can 
conclude that a xy = 0, implying that there is no Hall conductivity. 



The boundary conditions we impose on (30) are as follows. Near the black hole horizon, 
A x , (f ty and (fl y have the behavior 

(1 - z) ±l " /3 . (31) 

Requiring that the near-horizon modes of the gauge field and spin two field are falling into the 



horizon implies that we must choose the minus sign in (31). Equivalently, when computing 



the retarded Green's function for the charge current we should impose causal boundary 



conditions on the equations of motion (30). Although with these boundary conditions, 



ip zy ~ if* ~ (1 — z) lu> ^' A x , the Lorentz invariant quantity ^p*n U f fll/ is finite at z = zu- 



In principle, given the boundary conditions (31) we could integrate the differential equa- 
tions for A x , tp ty and yf ty numerically out to the boundary, located at z — 0, and look for 
solutions where the leading non-normalizable term in the near boundary series expansion of 
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ifty and <p* vanishes. From the relation between the leading and subleading terms of A x , we 
can read off the conductivity in the ^-direction. In practice, since the response of tp ty and 
(pl y to the electric field is linear, as is the response of A x , (p ty and <^* y to the source terms for 
ifty and (pi , we do not need a shooting algorithm to determine a. Instead, we can directly 
read off the xx component of the conductivity from the 3x3 matrix relating the expectation 
values of the boundary fields to the sources. The numerical value of the conductivity for a 
A = 4 condensate can be found in figure [2] 

It turns out that the (linear) conductivity is isotropic, in spite of the presence of a <i-wave 
condensate. Indeed, to get the conductivity in the y direction, we look at the effect of a 



7r/2 rotation of the condensate in equations (30). Such a rotation flips the sign of ip. By 



inspection of the equations, this is equivalent to flipping the sign of A x . Since flipping the 
sign of A x is equivalent to flipping the sign of both the electric field (the A^ term) and the 
current (the A^ term), the conductivity will be unaffected by such a change of sign. We 
conclude that the conductivity is proportional to the identity matrix. 

We note that the isotropy of the conductivity is only true at linear order. One can check 



that the non-linear version of (30) does depend on the angle between the electric field and 



the condensate, so the non-linear response of the material to an electric field is anisotropic. 

5 Fermions 

In this section, we compute a fermion spectral function in a holographic superconducting 
background. Our calculation is motivated by angle resolved photoemission spectroscopy 
(ARPES) experiments for d-w&ve superconductors. ARPES experiments reveal a great deal 
of information about the various phases of the cuprates and other superconductors. In these 
experiments, high energy photons striking the surface of a sample liberate electrons. The 
intensity of the produced electrons is a direct measure of the spectral function of electrons 
in the sample. Through this virtually direct method of measuring the electron spectral 
function, ARPES experiments provide an arena where experimental and theoretical results 
can be compared. See for example jlHl HD S2] for reviews of ARPES experiments. 

5.1 The fermionic equations of motion 

Consider a fermionic operator Oq, in the boundary theory dual to a bulk spinor \1/ whose 
action is 

S 9 = - J d d+1 Xv ^U + S hnd . (32) 
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Figure 2: (Color online) The real part of the conductivity as a function of frequency for 
a d-wave condensate of conformal dimension A = 4. As the temperature is decreased one 
observes a spike in the conductivity indicative of a bound state. This spike is localized at 
smaller values of u as the temperature is lowered. A second spike in the conductivity appears 
to vanish as the temperature is decreased. 



The Dirac Lagrangian has the standard form 



Cv = im (r»D„ - m c ) v 



(33) 



and the sign in (32) has been chosen to ensure that the spectral function is positive. The 



covariant derivative in (33) is given by 



'// '■'!> ^»IX,AQ*- l 1C,^-pu 



Dn, = d, L + -U,, A(j r 



(34) 



with lo the spin connection and ^ a Dirac spinor of charge q^. The bulk Gamma matrices 
P satisfy the Clifford algebra {F X ,P / } = 2g» u and P 4 " = r^Fl Vielbein indices are 
underlined. The dimension of Oq, is A# and is related to the mass of the bulk fermion m^ 
through 



d , 
A$ = — ± m^L . 



(35) 



The sign in (35) depends on the quantization scheme of the boundary theory and must be 



compatible with the unitarity bound A^ > (d — l)/2 [2H]. Because our spacetime has a 
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boundary located at z = 0, we entertain the possibility of a boundary action Sb n d- We will 
discuss the form of Sbnd once we have analyzed the near boundary behavior of \l/. 



The action (33), while familiar and natural, won't generate interesting dynamics. The 
fermions in (33) are coupled only to the gauge field and the metric and are not directly 
affected by the spin of the condensing field. Indeed, given the similarity between our con- 
densing spin two field and the Abelian Higgs model of [TU] . the behavior of the fermion 
correlator in the condensed phase will be identical to that of the s-wave superconductor 
studied in gHHI]. 

Thinking in the language of effective field theory, to take into account the effect of the 
ci-wave condensate we should add an interaction term to the fermion Lagrangian of suitably 
low dimension. For definiteness, in what follows we will work with d = 3. If we consider only 
interaction terms of mass dimension smaller than six, then we find three possible interaction 



terms: 
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%^| 2 W( Cl + zc 2 r 5 )$ , 

^^(c 3 + c 4 r 5 )Hk, (36) 

rfifl v WV^D u ^ + h.c. , 

where \l/ c = CT-ty* and the charge conjugation matrix C is defined by CT^C -1 = — T^ 1 ". 
The first term contributes to an effective mass for the fermion, without any angular depen- 
dence related to the <i-wave condensate. Since we do not expect any interesting physics to 
arise from it we will not consider it further. The second and third terms are not mutually 
compatible: the second term requires that the charge of the fermion, q^, equal the charge 



of the condensate, q, while the third term requires 2q^ = q. The second term in (36) will 



produce an s-wave gap similar to the one studied in [43J in a slightly different context. 



In what follows we will focus on the last term in (36) since it has the most interesting 



effect on the fermion correlator. Thus, the bulk fermion Lagrangian we consider is given by 

U = W (T^D^ - m c ) ^ + rf^WT^D^ - rj$T»D u (</V* C ) • ( 37 ) 

Since the phase of rj can be changed by a redefinition of the fermion field, we will set r] > 



1 There are other terms of dimension smaller than six that one can write down in the action. These are: 
i|t^| 2 *(c+icr 5 )*, (p*W(c + cT 5 )^ + h.c, ip* 2 W(c+cT 5 )^ + h.c, </?* ^r 5 !^* + h.c. where all the coupling 



constants are denoted by c. However, since our ansatz ( 16 ) satisfies the constraints ip = ip^ = and since 



v]/cpMv[/ — o due to Fermi statistics, none of these contribute to the EOM's of the fermion. 



18 



in the rest of this discussion. The Dirac equation following from (37) is 



= (P 4 /^ - m c )^ + 2ir]^^D u m c + ir)(ff,T^ c 



(38) 



It will be convenient to rewrite the equations of motion in terms of a rescaled four-component 
spinor £ = (—g ■ g zz ) x '^ which we will further decompose into two two-component spinors 
C = (Cij C2)- Since the Dirac equation involves both \l/ and \I/ C then, in Fourier space, it will 
couple modes with wavevector (u, k) to modes with wavevector (—a;, —k). Thus, we look for 



solutions to (38) which are of the form 



C = e- 



-iuit+ik-x /-(uj,k) 



p u >*)( z )+e 



iujt—ik-x ;-(—uj,—k) 



(39) 



where the notation ^'^ implies a single Fourier mode. We make the following choice of 
(real) bulk gamma matrices: 



r^ 



-ia 2 
ia 2 



a 3 
a 3i 



-prr 



ai 

CTi 



p2/ 







-io 2 



ia 2 



(40) 



and use C = T-. The Pauli matrices are denoted a%. The boundary gamma matrices 7 m are 
defined by the action of the Lorentz generators iT— on the positive eigenspace of T-, We 
find that 7' = — ia 2 , 7^ = Oi, j y — — 03 and C = —7*. With these definitions, the Dirac 



equation (38) reduces to 



= D (l) CP kx) + 2r ] {g xx f' 2 k x [p^-yyffiQ 



= D (2) CP kx) + 2 V (g xx f 2 k x [pxx-wtTiC. 



(—u>,—k x )* ■ /-{—U,—k x )*~\ 



(41) 



where we have aligned the momentum vector k along the x axis. The solution for arbitrary 
k can be recovered by using (17). The differential operators D( a ) are given by 



D 



(a) 



^a 3 d z - m c + (-l) a (w + q c A t )^g Tt a 2 + ik x y/^ ' a x . 



(42) 



The equations for £(~ w <~ fca; ) and for £(^>M* can De obtained from (41 ) by an appropriate sign 
flip or complex conjugation. 

We are interested in the retarded Green's function for the boundary fermion operator 
Oq, in the condensed phase at non-zero temperature. The prescription we use to study this 
fermion correlation function closely follows that of (13 EHJ H7J 08] which is based on the 
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work of [49]. In what follows we will summarize this prescription leaving a more detailed 



review to AppendixLD] Consider a solution to the Dirac equation (41 ) with infalling boundary 



conditions at the black hole horizon. We denote the coefficients of a near boundary expansion 



of the fermion fields which solve (41 ) as S and R where 



0{z) 

(o"i5') Q/ 



-m^L 



+ 



R 

P{z) 



y m ( L 



(43) 



The index a = 1,2 is a boundary spinor index, and the expansion above applies to both 
£(w,fc) anc j £(-u>,-A;)_ jj ere g j s identified with the source for the operator Oq, while R = (Oy). 
Since the Dirac equation (38) is linear, R^ u > k ^ will be linearly related to S^ u,k ' and S^~ u '~ k ' c : 



r£$ = MjS^ + MjS<f"'- %)c . 



(44) 



The prescription to compute the retarded Green's function is then [45] 



G r (uj, k) = -iMrf ■ 



(45) 



Typically the spectral function is defined to be the matrix quantity p(u), k) = (Gr(u), k) — 
G}j(a;, k))/2i. With a little abuse of terminology, we define the spectral function as 



p{uj, k) = Tr Im Gr(u, k) 



(46) 



where the trace is over spinor indices. 



5.2 The spectral function for fermion operators 



We obtained the spectral function (46) numerically using (45). After obtaining a numerical 



solution for the condensate (described in section 3.2) we solve the Dirac equation (41) and 



numerically extract the near boundary components S and R defined in (43). In what follows 



we focused on m^ = 0, leaving the exploration of the mass parameter to future work. 

Consider first a configuration in which the coupling of the fermions to the spin two field 
vanishes, rj = 0. As discussed earlier, since the fermions are not coupled to the condensate, 
the resulting Dirac equation is identical to the one studied in [4*3114*4"]. For configurations in 
which the near horizon geometry is asymptotically AdS and at finite q and zero temperature, 
refs. [HH HI] found that the spectral function for the fermions has continuous support inside 
the lightcone, \u\ > cm\k\, where cm is the speed of light in the infrared. For \u\ < cm\k\ 
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Figure 3: (Color online) The fermion spectral function for decoupled (77 = 0) massless 
fermions at T = 0.66T C . Red implies a large value of the spectral function and blue represents 
a value closer to zero. The red bands in the spectral function are associated with quasi-normal 
modes of £1 and (2. As described in the text, if the dispersion relation of (1 is 00 = E(k x ) 
then that of (2 is 00 = E(—k x ). 



the spectral function has delta-function support along codimension-one surfaces, which are 
in one-to-one correspondence with zero-modes of (. At T = one identifies the momentum 
\k\ = kp where a mode crosses the u = axis with a Fermi momentum in the bound- 
ary theory. As expected for an s-wave configuration, the spectral function is rotationally 
symmetric. 

As opposed to (121 HI] , we work at non-zero temperature and in the probe limit where the 
geometry is unaffected by the matter fields. As we approach T = the geometry approaches 
empty AdS and, in particular, its speed of light in the infrared is one|]_j We conclude that in 
the probe limit, and at low but non-zero temperature, there should be an approximate IR 
lightcone, given by \u\ > \k\, in which the spectral function has continuous support. Also, 
the zero-modes which will be observed at zero temperature should broaden into quasi-normal 
modes. These broadened quasi-normal modes can be seen in figure [3j 

Still working with rj = we focus, without loss of generality, on the k y = plane. 
The equations of motion for d and (2 decouple and we can separately discuss the quasi- 



11 We once again emphasize that when m 2 > the analysis carried out in |36[l37j implies that one can not 
consistently take q — > 00 and then T — > 0. Therefore, in this work we consider only T > 0. 
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normal modes of £i and £2- These quasi-normal modes appear as peaks in the components 
GR t n(u>,k x ) and G\r,22(w, k x ) respectively. Consider a single quasi-normal mode of, say, (j. 
We define the dispersion relation of this quasi-normal mode as the location of the maxima of 
G\r,ii and denote it by w = Ei{k x ). Generically there will be multiple quasi- normal modes 
so we add an index a to the dispersion relation to = E{" (k x ). Similarly, ( 2 will have quasi- 
normal modes with dispersion relation u> = E\ ? (k x )- For 77 = the equation of motion for 
(2 can be obtained from the equation of motion for (j by making the substitution (1 — > (T3C2 
and k x — > —k x . Therefore, for each quasi-normal mode of £i there is a corresponding quasi- 
normal mode of £2 and vice versa, such that E^ \k x ) = E[ a (— k x ). See figure 3 We note in 
passing that the quasi-normal modes of Q have dispersion relation —to = E\ ■ (— k x ). 

A non-vanishing, but small, rj/L in (41) couples the quasi-normal modes of ^^ and 



Q(-w-k x )* w j^] 1 an an gi e dependent strength, 9 being the angle between the condensate and 
the momentum vector. We define 9 = as the angle for which cp xy 7^ 0, tp xx - yy = and 
k = (k x , 0). The spectral function for this configuration is depicted in the left panel of figure 
4j At this angle Q u ' is coupled to Q w ' . Once the coupled modes intersect, eigenvalue 



repulsion will replace the intersection by a gap. This gap-generating feature of (37) is similar 
to the mechanism described in [43], where a gap was generated in an s-wave superconductor 
by adding a Majorana-like coupling to the fermion action. The set of lifted degeneracies 
includes all those located at u — 0. 

Next, consider a configuration at 9 = it/4 where (f xx ~ yy 7^ 0, (p xy = and k = (k x , 0) 
(this configuration is equivalent to one with k x = k y and tp xx - yy = 0). The spectral function 
for this configuration appears on the right panel of figure 4 We find that Q u ' is coupled to 
Q w ' and Q W ' " is coupled to Q "' ■ The set of lifted degeneracies is distinct from 

the one at 9 = 0. We call 9 = the anti-nodal direction and 9 = tt/A the nodal direction. 
A node is a point on the u = axis where two quasi-normal modes intersect. At 9 = 7r/4, 
Ci is not coupled to Q UJ ' and since only Q^' and Q UJ ' have dispersion curves 

that intersect along the u = axis, there will be no associated gap at such an angle. At 
intermediate angles 9 7^ 0, tc/A all the degeneracies are lifted. 12 This lifting is depicted in 
the central panel of figure [4| 

We note that once r]/L becomes too large our perturbation-theory based analysis breaks 
down and we find that a gap is generated at all angles. The critical value of 77 for which this 
happens depends on T/q^/x, A and m^L. 



12 To be more precise, we should write 9 ^ Q + n \ , 9 7^ f +n \ , because the spectral function is invariant 
under a n/2 rotation: The spin two field (p^ changes sign under a rotation of 7r/2, but this sign change is 
equivalent to changing the sign of r\ which in turn can be absorbed in a redefinition of the fermion field. 
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Figure 4: (Color online) The spectral function for the fermions evaluated for r] = 0.5 and 
T = 0.66T C . The red bands corresponds to a large value of the spectral function and the 
blue regions to values which are closer to zero. The left panel corresponds to the anti-node, 
the central panel to a 22.5 degree angle from the node and the value of the spectral function 
at the node is shown in the right panel. 





Figure 5: The Feynman diagrams associated with the interaction terms between the fermions 
and the spin two field in A. (a) rj*<p*^WV^D v ^; (b) rj^T^D" (if^ c ). 



There are other interesting features that come out of the analysis of the spectral function 
and that, surprisingly, compare well with properties of <i-wave superconductors. For instance 
the precise angular dependence of the gap is well fit by | cosfc x — cos/c y |, the nodes broaden 
into Fermi arcs in a particular temperature range, and Dirac cones are observed with Fermi 
velocities whose ratio can be tuned by dialing 1]. This analysis is performed in the companion 
letter [TTj. 



5.3 The physics of bulk fermions 

The bulk fermions are charged under a U(l) field and we can think of them as electrons 
and positrons. In the presence of a spin two condensate, the vertex tp* ^/ c T fl D u ^/ and its 
Hermitian conjugate lead to electron positron mixing. See figure [5j These interaction vertices 
have an interesting spin structure which we now elaborate. 

Consider the Dirac equation, D( a )( a = 0, in a limit (and region of space-time) where 
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k x = A t = and where the z- derivative can be ignored. In this limit, we can find approximate 
solutions to the Dirac equation of the form 

cr> « cr ,0) * - (_' .) , ct 0) - d" 0) * - ( ; ) , (47) 



where ui = m^/y/—g u . The Lorentz generator of rotations in the xz-plane is given by 




Thus, Ci ' has spin-up in the y-direction while Q^' has spin down in the y direction. 
In other words, in an approximate rest frame for the fermions, the quantum field Q u ' 
annihilates bulk spin up electrons in the y direction (for u > 0) or creates bulk spin down 
positrons (for ui < 0), while Q^' annihilates bulk spin down electrons or creates bulk spin 
up positrons. 

To gain further physical insight into the existence of these quasi-normal modes, consider 
a semiclassical picture of an electron with positive electric charge q^ > 0, positive spin in the 
y direction, and k x > 0. Such an electron is represented by Q w ' . This electron feels both a 
gravitational attraction to the black hole and an electrostatic repulsion from the electric field 
E z < 0. If the electric field becomes strong enough, a bound state appears in the potential 
for C{"' M 

If we boost to the rest frame of the electron, we observe a magnetic field B y < in addition 
to the electric field E z . This magnetic field will act to lower the energy of the electron because 
of the coupling with its spin. In other words, there is an additional shift to the energy level 
E\{k x ) of the bound state because of spin-orbit coupling. For the corresponding electron 
with spin down (Q^ ' ) the energy E2(k x ) is increased. If the potential, before considering 
spin-orbit effects, is not deep enough, the increase will lead to the loss of the bound state. 
However, as the electric field is increased, the potential well supports more and more bound 
states, and we see more and more alternating bands of d and ( 2 quasinormal modes. For 
large mass fermions, we expect the spin-orbit coupling to be a small effect, smaller than the 
spacing between the bound states in the potential, ensuring that the bands are alternating. In 
the small mass limit, our non-relativistic quantum mechanics intuition fails, but the numeric 
calculations of the previous section indicate that the bands continue to be alternating. 

The Fermi sea in the bulk consists of a series of concentric disks in the xy-plane in /c-space, 
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one disk per band of normal modes. Each disk possesses a Fermi surface at u = 0. In the 
limit in which k x is small and the rest frame calculation of the spin is a good approximation, 
as we move around one of these Fermi surfaces in fc-space, the spin of the electron rotates in 
the xy plane, always pointing tangent to the surface. In particular, for k x < and k y = 0, 
the role of £i and ( 2 are switched because the sign of the magnetic field is reversed if we 
boost in the opposite direction. 

From a bulk point of view, the coupling between the spin two field and the fermions 
produces new quasiparticles in the presence of a spin two condensate. In £;-space, at 9 = 0, 
electrons can mix with positrons of the same spin on opposite sides of the Fermi sea — (j with 
Q or £2 with (*. However, at 9 = 71/ '4 electrons mix with positrons with opposite spin on 
opposite sides of the bulk Fermi sea — Ci with Q or ( 2 with Q. This type of new quasiparticle 
with indefinite charge also appears in BCS theory. 

The spin structure of the bulk Fermi surface is crucial for the angle dependence of the 
gap. Naively, in the absence of a magnetic field, one would expect fermions with different 
spin to have the same energy at the Fermi surface. Because of the spin-orbit coupling, we 
find in fact that the energy levels are split. The dispersion curves for d and £* do not cross 
at the Fermi surface while the curves for (1 and Q do. Thus we find a node at 9 = ir/4 and 
a gap at 9 = 0. 

Do the fermions back react? Adding massless fermions to the system with charge q^, 
we would like to estimate semi-classically the total amount of charge contained in the Fermi 
sea. We restrict to d = 3 for this estimate, and we work locally, assuming the curvatures are 



small and space looks locally flat. With a potential of the form (21), the Fermi energy will 
be 



z / _ z 



E F = q c A tV /-g* = J-7jMc ( 1 " - J • (49) 

For a free massless gas of fermions, Ep = (37r 2 p) 1 / 3 where p is the local number density. 

Thus we have 

1 r / \n3 

(50) 



67T Z 



Z / 

- r fiq c 1 



The number density in the field theory from the fermions can be estimated by performing 
the integral 

n= r^/7 dz={v ^ 7T - 3) i^ zh - (5i) 
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The field theory charge density from A t on the other hand is fi/zh- To be able to ignore the 
back reaction of the fermions on the gauge field, we need 

tfzlql « 1 . (52) 

For this condition to be met at temperatures of order the phase transition temperature, we 
need A 2 q^/q 2 <C 1. We have chosen q^ = q/2 in previous sections, and so qA <C 1. 

Earlier, we mentioned that the probe limit involves taking q large and one might worry 
that the limit is incompatible with ignoring the back reaction of the fermions on the gauge 
field. In fact, we should be more careful. Although we don't have an action for the spin 
two field where we can rigorously take back reaction into account, we can look at the back 
reaction of the gauge field on the metric for temperatures close to T c . In this case, we know 
the fully back reacted solution involves a Reissner-Nordstrom black hole in AdS. Placing a 
1/2k 2 in front of the gravitational action, we find that the back reaction modifies the warp 
factor: 

f(z) = l-(j-) (l + Q)+(j-) Q, (53) 

where Q = z 2 .k 2 h 2 /2L 2 . In the normal phase the probe limit condition is thus 

KZhfi < L . (54) 



Close to T c , using (87), the inequality reduces to q ^> ran. Thus, if we can tune ran to be 
suitably small, then there is an intermediate coupling regime run <C q <C 1/mL where the 
probe limit is valid and fermionic back reaction can be ignored. Whether such a regime can 
be realized in a stringy construction is another question. 



6 Outlook 

In this work we have discussed the construction of a holographic <i-wave superconductor. 
As emphasized several times in the text, our construction is not ideal. There are well- 
known difficulties in writing down an action for a charged spin two field propagating in 
a curved spacetime. These difficulties lead to non-hyperbolic and non-causal behavior of 
the spin two field. In principle these problems could be circumvented when considering the 
phase diagram of the theory and in studying the fermion correlators in the condensed phase 
if one works with a very massive spin two field in the bulk (or, equivalently, a spin two 
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operator of large conformal dimension in the boundary). Numerically, it is challenging to 
construct backgrounds for fields of very large mass and it would be worthwhile to improve 
on our current results where we studied condensates of dimension less than or equal to six. 
Alternately, one might try to use analytic techniques, perhaps similar to the ones used in 
appendix [C] or in [HU] to construct these solutions. 

It would be, of course, much more satisfactory to have an action that produces hyperbolic 
and causal equations of motion with the correct number of degrees of freedom. Unfortunately, 
such an action is not presently known. There are, however, several different directions one 
could try to pursue. In [5T] a consistent action for a spin 3/2 particle propagating in a 
constant electro-magnetic field and flat spacetime was constructed. Similar techniques might 
be used to obtain an improved form of the Argyres-Nappi action [22] which would have the 
correct number of propagating degrees of freedom in dimensions other than 26. A different 
approach that would lead to a causal action for a spin two particle would be to consider a 
Kaluza-Klein compactification of Einstein gravity on a manifold with a U(l) isometry. The 
drawback of such a method is that, generically, such a compactification would lead to an 
action describing a whole tower of massive spin two fields. A final method which might prove 
useful would be an expansion of our proposed action in powers of F^ v and its derivatives, 
requiring that the resulting equations of motion are causal and have the right number of 
degrees of freedom at each order of the expansion. 

Even without a full-fledged causal action for the spin two field, there are certain properties 
of the condensate that may be probed. In this paper we discussed only solutions where tp xy 
and (fxx-yy are real. By introducing a relative complex phase between these modes, one 
obtains an exotic d + id order parameter that breaks time reversal invariance. This d + id 
superconductor will probably have a non-vanishing Hall conductivity (with spontaneous 
currents along the boundaries of the sample) in the absence of a magnetic field. Similar 
phenomena were reported in a holographic context for p-wave superconductors in [52, 153] . 

In addition to the conductivity, we studied the effect of the spin two condensate on probe 
fermions. We saw how the angle dependence of the gap emerged from simple qualitative 
arguments involving degenerate perturbation theory and spin-orbit coupling. Given the dif- 
ficulty of solving the Dirac equation, it would be interesting to see if semi-classical techniques 
such as the WKB approximation could be used to further elucidate the physics of fermions 
in these black hole backgrounds. 
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A The spin two Lagrangian 

It is difficult to write down an action for a massive charged spin two field in a curved space- 
time. Most of the difficulties are associated with irregularity of the spin two Lagrangian 
which leads to a constrained system of equations. In what follows we will begin by review- 
ing the Fierz-Pauli Lagrangian for a neutral, non-interacting spin two field in Minkowski 
space [TB]. Many of the strategies used in this simple case can be abstracted to the more 
complicated example of interest. 

A.l Neutral field in flat spacetime 

The EOM's that follow from the Fierz-Pauli Lagrangian (pi) are 

E w = VV,a ~ ^m," _ V^ + <£V - ™ 2 </V _ V^V^a _ </>A, A _ m V) = • (55) 

Inspection reveals that the d + 1 equations E^ t = do not involve second order time deriva- 
tives. Thus, the Lagrangian ([I]) is irregular and the equations of motion E^ t = are actually 
constraintsp] The existence of constraints is expected since the symmetric tensor ip^ u naively 
has d + 2 more degrees of freedom than a massive spin two particle. In order to fix the values 
of the extra d + 2 fields and their momenta we need 2(d + 2) constraint equations. The 
divergence of E^ u involves only first derivatives and yields d + 1 further constraints: 

E^ = m 2 (<p >v -(p„)=0. (56) 
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For a brief explanation and classification of constrained systems see e.g. [53] . 



A penultimate constraint can be obtained by using the trace of the EOM's, 

E£ = (d- l)(<p x x - tp x ) + dm\ = , (57) 

together with E pu ^" = which yields 

<p = 0. (58) 

The last constraint is ip y t = 0. In this simple case, the EOM's and constraints can be 
rewritten in the form (J2J), which is the form of Fierz. 

A. 2 Charged field in curved spacetime 

Having tackled the Fierz-Pauli action, we are better prepared to deal with the more involved 
case of a charged massive spin two field in a curved background. Our treatment follows 

pi using. 

Let us write down the most general d + 1 dimensional action, quadratic in the spin two 
field and of scaling dimension d + 1: 

C = -\D p if^\ 2 + 2\D^\ 2 + \D^\ 2 - [D^ v D v ip + h.c] - m 2 (|w| 2 - M 2 ) 

+ c x Rp,pW*">(p vX + c 2 R^ X ifl + c 3 [e ie R^<p*^<p + h.c] + c 4 i?|</v| 2 + c 5 R\<p\ 2 

+ iceqF^cp^tfi . 

(59) 

The Q are assumed to be real. The terms in the first line of (59) are fixed by requiring 



that ( 59 ) reduces to (HI) when the metric is Minkowski and when the charge of the spin two 



field vanishes. The terms on the second line of (59) involve all possible curvature invariants. 
They vanish on a flat spacetime. There is only one possible invariant we can construct from 
F pv that is quadratic in the ip pv and of scaling dimension d+1. 

We point out two subtleties in this Lagrangian. First, we can rule out a contribution of 
the form 

i [D^ v D v ip -h.c] (60) 

because in the q — > limit such a term would mix the real and imaginary parts of (p^ instead 
of leading to a Lagrangian for two, free, real spin two fields. Second, there is a potential 
ambiguity that involves a combination of the first and second lines: the expression 

v^ ( - D p( p*^D x <f^ + D x v *^D pVx , - R^p^V + RxpV*^ + iF^^vl) (61) 
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is a total derivative]^ We have resolved this ambiguity by setting the {D\tp*^ v )(D tl Lp Xl/ ) term 
in the first line to zero and allowing for c\, C2, and cq in the second and third lines. 

We will fix the q and 9 by requiring the existence of 2(d + 2) constraint equations. We 
begin by computing the EOM's: E pu = 0. Just as in flat space, the equations E pt = do 
not involve second derivatives with respect to time and constitute d + 1 constraints. Taking 
the divergence of the EOM's, D p E pu = 0, we find d + 1 additional constraints. 

To find the two remaining constraints, we pursue the same strategy used in flat space: 
we would like to combine the trace equation E^ = with D fJ- D l/ E pu = in such a way as 
to eliminate the second derivative terms, leaving at most first derivatives of the fields. The 
trace of the EOM's, E£ = 0, reads 

[dm 2 +(c 3 + c 4 + (d+l)c 5 )R]^+(c 1 + c 2 + (d + l)c 3 )R^^ pu = (d-l)(n^-D^ p ) , (62) 

while D tl D u E pv = can be written in the form: 

-m 2 (D <p - D^) = iq(l + c 6 )F pa D p <p a + ( Cl - 2)R" oa,p D lt D v ip a0 + c 2 W a D^ a 

+ (1 + c 3 e ie )R» a D p D aV + c 3 e- ie R a P □ <p a0 + Ci RD p ^ p + c 5 R □ <p 
+ ... (63) 

where the ellipsis denotes terms that involve at most single derivatives of (p pv . 



Since there are seven independent terms on the right hand side of (63) and only six 
coefficients q, we cannot make the right hand side vanish and we need to invoke an additional 
assumption about the metric g pv . If we assume that the metric satisfies the Einstein condition 
(|6| then eq. (63) becomes 



m 2 (D if - D^) = iq(l + c 6 )F pa D pVa + ( Cl - 2)R^D fi D v <p a p 
2A 



+ 



d-1 



(c2 + (d+l)c 4 )D t *<p fl +(l + 2c 3 cos9 + (d+l)c 5 )nv + ■■■ (64) 



The Cj's can now be chosen so that the right hand side is equal to (m 2 , — m 2 ) (□ tp — D^ip^) 
with m,Q a real constant: 

ci = 2, c 6 = -l, -(c2 + (d+l)c 4 ) = (l+2c 3 cos0+(d+l)c 6 ) = — — (m 2 -m 2 ) . (65) 
As suggested by the appearance of m 2 , above, this one parameter family of solutions can be 



l Our convention for the Riemann tensor is that [V^, V y ]l/ A = R x ppv V p . 
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absorbed into a rescaling of the mass. We are led to the following Lagrangian density: 

C = - |ZV/Vl 2 + 2|iWf + \D^\ 2 - [D^ u D uV + c.c] - m 2 (|yv| 2 ~ M") 

1 (66) 

+ 2R llvpX ip*>»tp 1 ' x - R^W - j^RW? ~ tqF^Vx ■ 



With the addition of a Maxwell term, (66) becomes (J3J). 

Before we proceed with the last two constraints, it will be useful to write down the 



equations of motion, their trace, their first divergence and their second divergence after (65) 
have been implemented. The equations of motion are 



E^ = (□ -m 2 )yv - 2D {p <p u) + D ip D u) ip - g pu [(□ -m 2 )ip - D p <p p ] 



2R f j,pi,\ip p — g^u 



pX ~" ^^P-^(F, P ^ + F up ^)=0. 



(67) 



The trace equation (62) reduces to 



EH 



dm 



d-1 

d + 1 



R 



v-{d-i)(n(p-D^) =0. 



(68) 



The first divergence can be written as 



D P E 



\AV 



-m 2 ((pv - D v tp) +iq 



3 
2 



1 .. 3„ „„ 1 



- t F v (p + -F^D'ip 



{D p F pv )v 



w 



0, (69) 



where we have employed the notation F v = D^F pv . The second divergence can be written 

as 



D V D^E, 



flV 



-m 



(D»<p u -nv)+iq {D»F v »){D v (p pp ) - F^ p + F»D^ 



q 2 F^F up ^ + -q 2 F pv F^ ¥ , = 0. (70) 



Our penultimate constraint is obtained by combining the second divergence (70) with 



the trace equation (68): 



— m 



dm 2 



R 



d-1 d + ll 



(p = iq 



(W)(A,vv) 



- F^ p + F»D^ 
3 



q 



F^Fytft - - F, V F^ 



(71) 
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Unlike the q = case, the constraint (71) is not a purely algebraic constraint on ip. Instead, 



it involves first derivatives of (p^ multiplied by first derivatives of F^ v . Thus, it is not 
immediately clear that taking a time derivative of this constraint equation will yield the 



final constraint. Acting on (71) with D l we find 
dm 2 



m 



R 



D l (p = iq 



(D J FnD t DM + (D l F^D t DM-5M) + ... , (72) 



.d-1 d + ll 
where here and below the ellipsis denotes terms that do not involve second order time 



derivatives of <^ M „. We now use the equations of motion (67) and its divergence (69) to 
eliminate these second order time derivatives. The space-space components of the EOM's 
give 

= E) = &D t {$ - S^t) + . . . , (73) 



which we can use to eliminate the second term on the right hand side of (72). Taking the 
time derivative of the divergence of the EOM's we get: 



= D l D^ 



-m 



3iq 



D'Dw) + -* [F'jD'DM + FM>* A(?J - <*>£)] + • • • 



(74) 



Using again (73), we can eliminate the second term in the square brackets. Finally, for 
generic values of F^ we can invert the matrix 



m 



3 2m 2 j 



(75) 



and solve (74) for D D t (p- in terms of first order time derivatives of (p^ u . With this solution 



in hand, we can eliminate also the first term on the right hand side of (72), and reduce it 



to the last constraint equation. As pointed out in [56] , there will be choices of Fij for which 
this matrix is not invertible, so that we fail to find a constraint. However for qF l j/m 2 <C 1, 
which is one of the requirements in (J7l), the matrix is invertible. 

Given the possible non-existence of the last constraint, one wonders if our Lagrangian may 
have other pathologies. Indeed [571 EI] discovered that for generic values of F^, the EOM's 
are either non-hyperbolic or non-causal. Failure of hyperbolicity (which is the condition 
to have a well-posed Cauchy problem) is associated with electric fields, and it is avoided 
by taking F l j smaller than a certain bound. On the other hand, non-causal behavior is 
associated with both electric and magnetic fields, and appears for arbitrarily small values of 
such fields. Fortunately this effect appears at an order specified by (m) and is small when 
the field strength is small. In a frame where F l j is small this pathology can presumably be 
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corrected by adding terms to the Lagrangian that are higher order in F^. Such an expansion 
was carried out in a different context in [58] . 

B Stability bound for higher spin fields in AdS 

Breitenlohner and Freedman found that there is a lower bound on the mass squared of a free 
scalar field propagating in AdS^+i in order for it to be stable and non-tachyonic 



m 2 L 2 > -^ . (76) 

- 4 

Their argument is based on the requirement that the energy functional for the scalar be 
positive definite. For fermionic fields, reality of the action requires m to be real so that 
m 2 > 0. We find that, for bosonic fields of spin one and higher, the stability bound on the 
mass squared is also expected to be 

m 2 > . (77) 



In what follows, we illustrate ( 77 ) for the simple case of a massive vector field. 



Consider a massive vector field, whose action is given by 

S = \j d d+1 x v^ ( - i V" - m 2 A^ , (78) 

propagating in global AdS^+i with line element 

ds 2 = (-dt 2 + dp 2 + sinh 2 p dtt^/cosh 2 p . (79) 

The energy functional 



E= I 'd d x^T ^ 
for the spin one field can be constructed by integrating the bulk stress-energy tensor 



JO) 



T --h^ (81) 

over a spacelike slice orthogonal to the timelike Killing vector ^ = (1, 0, . . . , 0). Distinguish- 
ing the time component (p = 0) from the spatial ones (p, v = i,j), one finds that the energy 
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functional (80) takes the form 



E = - d d x tanh^ 1 p \g ij F oi F 0j + m 2 A 2 + — F^F* + -^— AA 1 } . (82) 

2 J I 2 cos 2 p cos 2 p ) 

For m 2 > the energy functional is manifestly positive definite. For m 2 = the energy 
functional and action are those of a massless vector field, and are both invariant under gauge 
transformations. Solutions with F^ v = have zero energy, and so the energy functional is 
semi-positive definite. Interpreting the energy as the norm in the Hilbert space of states, 
one finds that such states are null and therefore not physical. When m 2 < solutions with 
Ffj, u = can have negative energy. The kinetic energy of these solutions is zero and the 
(negative) mass term in the energy functional is responsible for the energy being negative r°\ 

The crucial point in the above computation is that the action is the sum of two terms: 
a kinetic term that is invariant under gauge transformations, and a mass term that is not 
invariant under gauge transformations and gives rise to a term of definite signature in the 
energy functional. Shifting our attention to higher spin bosonic fields, whenever the action is 
the sum of two terms — a kinetic term and a mass term — with the properties just described, 
stability of the energy functional should imply m 2 > 0. Note that bosonic fields of spin 
larger than two require auxiliary fields; therefore our argument regarding positivity of the 
energy functional requires improvement. 

Positivity of the energy functional is meaningful only when it is finite and can be used 
to bound A. Considering the neutral spin two field with BGP Lagrangian [19] (J3J), with 
the addition of suitable boundary terms, the energy functional (or, equivalently, the on-shell 
action) converges if ip^ vanishes at large p faster than (coshp)^~. Since (p mn ~ cosh ~ 2 p, 
A should be positive for d > 3. Given the relation m? = A(A — d) and the bound m 2 > 0, we 
conclude that in fact A > d. As opposed to the scalar case where the leading and subleading 
terms of a near-boundary series expansion can change roles [M] , the spin two field has only 
one possible quantization scheme. 



15 One might be worried that the constraint equation D^A^ = prohibits solutions of the form F^ = 0. 
However, since (82 1 does not depend on OqAo nor d p A p , at a given time we can construct modes whose only 



non-trivial component is A p (p), vanishing at large p arbitrarily fast, and with doA p = and OqAo fixed by 
the constraint. 
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C Approximating T c in the classical limit 



In section 3J2 we have argued that the system of equations governing the phase transition for 
a massive charged spin two field is identical to that for a massive charged scalar field. When 
the mass of the scalar field becomes very large we can treat the system semi-classically 

Consider a massive, charged particle in a charged AdS black hole background. Such a 
particle experiences an attractive gravitational force toward the black hole and a repulsive 
electrostatic force towards the asymptotically AdS boundary. For a high temperature black 
hole the gravitational force will always dominate over the electrostatic repulsion, and there 
will be no equilibrium configuration. However, as the black hole cools, the gravitational force 
close to the horizon becomes weaker, and one finds that the particle will settle in a global 
minimum of its potential, away from the black hole horizon. We identify the temperature 
at which this global minimum appears with the temperature of the phase transition for the 
scalar condensate. 

The action for a massive, charged point particle is given by 



S PP = —m \ dr — q A^dx^ . (83) 



The effective potential for a static configuration in the background (21 ) and in d = 3 bound- 
ary dimensions is given by 



mL y/1 - C 3 | C ~ 1 , 

z h C z h 



V pp(0 = — , + ^— ft (84) 



where we have defined £ = z/zh and fi = q^iZh- We obtained (84) by evaluating the point 



particle action (83) on a solution where the point particle is stationary. The minimum of 



the potential is located at a C, m in for which 

V^(Cmin) = • (85) 

Since V(l) = this minimum will be a global minimum only when p, is large enough so that 

V pp {( m in) < . (86) 
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Figure 6: (Color online) The critical temperature for generating a holographic <i-wave or 
s-wave condensate as a function of the dimension of the condensing operator. The numerical 



value for the critical temperature is plotted in solid blue. The approximation (87) is plotted 
as a dashed black line. 



A quick computations shows that (86) and (85) are satisfied whenever \x > \x c where 



fjL c 



+ — (152072-1224^) + -(19009 + 153^) mi~ 3.68 mL . (87) 
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In figure 6 we compare the expression (87) with a numerical evaluation of the normal modes 



of a condensing scalar field. Note that the numerical value for T c is slightly higher than the 



classical approximation (87). This is, perhaps, expected: in the quantum theory the ground 



state wavefunction has some kinetic energy so the minimum of the potential ought to be 
slightly lower than zero. 



D Fermion Green's functions 

We define the retarded Green's function for fermions in analogy to the Green's function for 
a scalar operator: 

G R (t,x) =iQ{t)({O 9 {t,Z),O\ f {0)}) . (88) 

In order to derive the retarded fermionic Green's function from AdS/CFT, we must first 
obtain the explicit form of the boundary terms in the fermion action Stmd- The variation of 



the Lagrangian (37) which leads to the Dirac equation (38) also results in a boundary term 
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of the form 



5S } 



bulk 



-/ / d 3 x(R5S-S5R) 



(89) 



where R and S were defined in (43) and are three dimensional spinors. In minimizing 



the action, we would like to vary Sy such that 5S vanishes on the boundary while SR is 
unconstrained. Thus, we need to add to the action a boundary counter term in order to 
have a well defined variational principle. An appropriate boundary action is 



& 



bnd 



i \ d 3 xSR = -- 



J d 3 x^^ (vv + ttn„r"tt) 



(90) 



where n^ is a unit vector normal to the boundary, pointing toward positive z. The total 
variation of the action is then 



5Sy = <5Sb U ik + £Sbnd = — if d 3 x (8S R + R 5S) 



(91) 



which is independent of 5R as required. 

Having found the appropriate boundary counter term, we return to a construction of the 



fermionic Green's function. Instead of thinking of (91 ) as a result of varying the action, we 



will now think of it as a perturbation to the action and ask how the fermionic operators 
respond at linear order to the presence of a 5 S ^ 0. From the theory of linear response, an 
expression that describes the change in the expectation value of an operator O due to the 
addition of a term 5H to the Hamiltonian is 



6(0(x)) = -i \ dt'e(t-t')([0(x),SH(if)]) . 



(92) 



For the fermionic operator R, we have J dtSH = —55*^ and 



S(R(x)) 



d 3 x'Q(t - t') [{{R(x), R(x')}}5S(x') + ({R(x), R c (x')})5S c (x')] (93) 



-i / d a x 



3 ' 



G R {x - x')Y5S{x') + G c R {x - x')Y5S c {x') 



(94) 



The Fourier transform of this relation is the multiplicative identity 



R (u,k) = _i(G fl (jfc) 7 * 1 s (w, * ) + G c R (k) 1 t S { - UJ '^ k)c ) 



(95) 



where Gn{k) is the Fourier transform of Gr(x). Comparing (95) with the bulk computation 



(44), we deduce the form of the retarded Green's function (45), consistent with the result of 
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Before wrapping up this section, we would like to check that we have chosen the sign 
of the fermionic action such that the spectral function p(k) = Trim Gn{k) is positive. For 
k = (cu, 0) and m^ = the dimension six operator in (37) vanishes and the Dirac equation 
can be solved exactly. We find that Ai(u, 0) = 7* and therefore p(cu, 0) = 2. 



This result for p(co,0) suggests the following sum rule 



/ duo (p(u, k) - 2) = 



(96) 



which we have been able to verify numerically for a few generic values of k. Analytically, we 
have found that for r\ = 0, p(u), k) = 2 + 0(l/u 2 ), although we believe the result holds more 
generally. Moreover, we expect that Gr is analytic in the upper half of the complex u plane. 
Thus one ought to be able to demonstrate the sum rule through a contour integration. 
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